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Motivated by the recent results from the PAMELA and ATIC, we study the 
cosmic-ray electron and positron produced by the decay of gravitino dark mat- 
ter. We calculate the cosmic-ray electron and positron fluxes and discuss im- 
plications to the PAMELA and ATIC data. In this paper, we will show that 
the observed anomalous fluxes by the PAMELA and ATIC can be explained 
in such a scenario. We will also discuss the synchrotron radiation flux from the 
Galactic center in such a scenario. 
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1. Introduction 

In astrophysics, the existence of dark matter (DM) is almost conclusive. Ac- 
cording to the recent survey of WMAP, 1 it accounts for 23 % of the total 
energy density in the universe. In the standard model of particle physics, 
however, there does not exist candidate for DM, which is one of the reasons 
to call for beyond the standard model. Supersymmetry (SUSY) is a promis- 
ing model which can give an answer to the question; in the framework of 
SUSY, lightest superparticle (LSP) is a viable candidate for DM. 

The fluxes of high energy cosmic rays give information about the prop- 
erties of DM. In the recent years, accuracy of the measurements of the 
fluxes have been significantly improved. In particular, anomalous signals 
are reported by PAMELA 2 and ATIC 3 in the observations of cosmic-ray 
e ± . The PAMELA and ATIC results have attracted many attentions be- 
cause the anomalies may indicate an unconventional nature of DM. In fact, 
there have been a sizable number of DM models are proposed to explain 
the anomalies after the announcements of the PAMELA and ATIC results. 
Among them, especially in decaying DM scenarios, the observed anomalies 
can be well explained with the appropriate choice of the lifetime of DM 
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especially in leptonically decaying scenarios. (For works calculating cosmic- 
ray e ± , see 4,5 and references therein.) 

In usual supersymmetric scenario, i?-parity conservation is assumed, 
which protects LSP from decaying into standard model particles and makes 
it a viable candidate of DM. If we consider the case that i?-parity is violated, 
LSP is no loner stable; however, if i?-parity violation (RPV) is weak enough, 
the lifetime of the LSP can be much longer than the present age of the 
universe and LSP can play the role of DM. 6 In addition, when the order of 
the RPV is properly chosen to give the lifetime of O(10 26 sec), produced 
cosmic-ray positron gives excellent agreement with PAMELA data. 4 

On the other hand, synchrotron radiation from the decay of DM may 
give constrains directly to scenarios explaining the PAMELA and ATIC 
anomalies. Since DM decays into energetic e under the magnetic fields in 
our galaxy, synchrotron radiation is inevitably induced. Importantly, the 
WMAP collaboration has observed the radiation in the whole sky, so that 
the observation gives constraints on the scenarios of the e ± production due 
to the decay of DM in the Galactic halo. 

In this paper, we consider gravitino (donated as tp^) LSP in RPV. In the 
scenario, we calculate cosmic-ray and synchrotron radiation flux induced 
by them, paying particular attentions to PAMELA and ATIC anomalies. 
We will see that the PAMELA and ATIC anomalies are simultaneously 
explained if the lifetime of the gravitino DM is O(10 26 sec) and the mass is 
^1 — 2 TeV. 4,5 In addition, synchrotron radiation from the Galactic center 
is comparable with or smaller than the observation. 7 

2. The Scenario and Model Framework 

In this section, we briefly explain the cosmological aspects and the model 
framework of ^-DM scenario in RPV. With RPV, the t/v LSP becomes 
unstable and energetic positron can be produced by the decay. Even if 
the W is unstable, it can be DM if the RPV is weak enough so that the 
lifetime of the gravitino t 3 / 2 is much longer than the present age of the 
universe. 6,8 In fact, such a scenario has several advantages. In the W-LSP 
scenario with RPV, the thermal leptogenesis 9 becomes possible without 
conflicting the big-bang nucleosynthesis constraints. In addition, the fluxes 
of the positron and 7-ray can be as large as the observed values, and the 
anomalies in those fluxes observed by the HEAT 10 and the EGRET 11 ex- 
periments, respectively, can be simultaneously explained in such a scenario 
ifT 3/2 ~O(10 26 sec). 12 - 13 

Here, let us consider the bi-linear RPV interactions. Using the bases 
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where the mixing terms between the up-type Higgs and the lepton doublets 
are eliminated from the superpotential, the relevant RPV interactions are 
given by 

£rpv = BiLiH u + m\ iHd LiH2 + h.c., (1) 

where Li is left-handed slcpton doublet in i-th generation, while H u and 
Hd are up- and down-type Higgs boson doublets, respectively. Then, the 
ipft decays as i/v — > ifW^, ViZ, Vih, and 1^7, where if and Vi arc the 
charged lepton and the neutrino in i-th generation, respectively. Taking 
account of all the relevant Feynman diagrams, we calculate the branch- 
ing ratios of these processes. 12 When the gravitino mass to 3 / 2 is larger 
than mw, the dominant decay mode is W ~~ * lfW T . In such a case, 
we see t 3 / 2 ~ 6 x 10 25 sec x (Ki/l0~ m )~ 2 (m 3 / 2 /l TcV)~ 3 , where Ki — 
(Bi sin/3 + m 2 f cos/3)/to? is the ratio of the vacuum expectation value 
of the sneutrino field to that of the Higgs boson, with tan/3 = (H^)/(H^), 
and rriOi being the sneutrino mass. Thus, r 3 / 2 is a free parameter and can 
be much longer than the present age of the universe if the RPV parameters 
Bi and m ~ are small enough. 



3. Electron and Positron Fluxes 

Let us first summarize our procedure to calculate the fluxes $ e ± . (For 
detail, see. 4 ' 5 ' 12 ) We solve the diffusion equation to take account of the 
effects of the propagation of e ± . The energy spectrum of the e ± from DM 
f e ± (E, r) evolves as 14 

^ = K {E)W 2 f e± + A [b{E)f e± ] + Q. (2) 

The function K is expressed as K = KoEq cV , 15 where EceV is the energy in 
units of GeV, while b = 1.0 x 10~ 16 xEq eV GeV/sec. In our numerical calcu- 
lation, we use the following three sets of the model parameters, called MED, 
Ml, and M2 models, which are defined as (5, if [kpc 2 /Myr], L[kpc]) = 
(0.70, 0.0112, 4) (MED), (0.46, 0.0765, 15) (Ml), and (0.55, 0.00595, 1) (M2), 
with R — 20 kpc for all models. Here, L and R are the half-height and the 
radius of the diffusion zone, respectively. The MED model is the best-fit to 
the boron-to-carbon ratio analysis, while the maximal and minimal positron 
fractions for E > 10 GeV are expected to be estimated with Ml and M2 
models, respectively We found that the MED and Ml models give similar 
positron fraction, so only the results with the MED and M2 models are 
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where tdm is the lifetime of DM. In the above expressions, [dN e + /dE]d cc 
is the energy distributions of the e ± from single decay processes, respec- 
tively, and are calculated by using PYTHIA package 16 for each DM can- 
didate. In addition, pdm is the DM mass density for which we adopt the 
Navarro-Frank- White (NFW) mass density profile: 17 Pnfw(x) = PQ r Q(r c + 
r ) 2 /r(r c + r) 2 , where p Q ~ 0.30 GeV/cm 3 is the local halo density 
around the solar system, r c ~ 20 kpc is the core radius of the DM pro- 
file, tq ~ 8.5 kpc is the distance between the Galactic center and the solar 
system, and r is the distance from the Galactic center. 

Once / e ± are given by solving the above equation, the fluxes can be ob- 
tained as [<£> e ± (E)] DM = 4^r/ e ± (E, Xq), where Xq is the location of the solar 
system, and c is the speed of light. In order to calculate the total fluxes 
of e ± , we also have to estimate the background fluxes. In our study, we 
adopt the following fluxes for cosmic-ray e produced by collisions between 
primary protons and interstellar medium in our galaxy: 14 ' 18 [$ e -]BG = 
0.1625^7(1 + ll£ G c 9 v + 3.225&#) + 0.70£° e 7 v /(l + 110^ + 600^ + 
580-Eq 2 v ) GeV 1 cm -2 sec -1 str -1 for the electron, and [$ e +]BG = 
4.5^ C 7 V /(1 + 650^ C 3 V + 1500££ 2 V ) GeV" 1 cm" 2 sec" 1 str" 1 for the 
positron. 



4. Synchrotron Radiation: formalism and the observation 

In this section, we first show the formalism for calculation of synchrotron 
radiation flux. (The detail is in. 7 ) Then, we address the implication of 
the present observation of synchrotron radiation from the Galactic center 
region. 

Synchrotron radiation energy density per unit time and unit frequency 
is expressed as 

L v {x)= I ' dEV(v,E)f e ±(E,x). (4) 

Here, Viv, E) is synchrotron radiation energy per unit time and unit fre- 
quency from single e^- with energy E. Adopting the Galactic magnetic flux 
density of B ~ 3 /xG, we can see that the synchrotron radiation in the 
observed frequency band of the WMAP (i.e, 22 — 93 GHz) is from the 
with the energy of E ~ 10 — 100 GcV. For the e ± in such an energy range, 
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f e ± can be well approximated by 




1 Pdm(x) Y e ±(> E) 

TDM WDM HE, x) 



(5) 



where F e ±(> E) = /~ dE'[dN e ± /dE'] dcc . Thus, we use f ( e ± > in the cal- 



culation of synchrotron radiation flux. Here, we note in this formula that 
we take into account the effects of both synchrotron radiation and inverse 
Compton scattering for energy loss rate as b(E,x) = P S ynch + -Pic- This 
is because, in the Galactic center region, the inverse Compton scattering 
in the infrared 7-ray from stars becomes the dominant energy-loss process, 
thus it can not be neglected. 

In order to calculate the observed radiation energy flux, we integrate 
L v {x) along the line of sight (l.o.s.), whose direction is parametrized by the 
parameters 9 and 4>, where 9 is the angle between the direction to the Galac- 
tic center and that of the line of sight, and <f> is the rotating angle around 
the direction to the Galactic center. (The Galactic plane corresponds to 
(f) = and it.) Then, the synchrotron radiation flux is given by 



Notice that, adopting the approximation of the constant magnetic flux in 
the Galaxy, the line of sight and energy integrals factorize. 

Radiation flux from Galactic center region has been observed by the 
WMAP for frequency bands of 22, 33, 41, 61, and 93 GHz. 19 - 20 Since then, 
intensive analysis has been performed to understand the origins of the ra- 
diation flux. (For recent studies, see. 19-21 ) Most of the radiation flux is 
expected to be from astrophysical origins, such as thermal dust, spinning 
dust, ionized gas, and synchrotron radiation, which have been studied by 
the use of other survey data. 22 With the three-year data, the WMAP col- 
laboration claimed that the flux intensity can be explained by the known 
astrophysical origins. 23 On the contrary, Refs. 19,21 also studied the WMAP 
three-year data, and claimed that there exists a remnant flux from un- 
known origin which might be non-astrophysical; the remnant flux is called 
the "WMAP Haze" . However, no clear indication of the WMAP Haze from 
unknown source was reported by the WMAP collaboration after five-year 



The existence of the WMAP Haze seems still controversial, and the 
detailed studies of the WMAP Haze using the data is beyond the scope of 
our study. Here, we adopt the flux of the WMAP Haze suggested in 21 (i.e. 
0(1 kJy/str)) as a reference value. 





data. 
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Fig. 1. Positron fractions for the case that i/> M dominantly decays to (a) the first- 
generation lcpton in MED model and (b) the second-generation lepton in M2 model. 
Here, we take m 3/2 = 300 GeV, 600 GeV, and 1.2 TeV (from left to right), with 
r 3/2 = 2.0 X 10 26 sec, 1.1 X 10 26 sec, and 8.6 X 10 25 sec (9.3 X 10 25 sec, 5.8 X 10 25 sec, 
and 5.0 X 10 25 sec) in (a) ((b)), respectively. Dot-dashed line is the fraction calculated 
only by the background fluxes. 

5. Numerical results 

First, we show the numerical results of the positron fraction. For simplic- 
ity, assuming a hierarchy among the RPV coupling constants, we consider 
the case where the i/v decays selectively into the lepton in one of three 
generations (plus W ± , Z, or h). In Fig. 1, we show the positron fraction 
for the case that the tp/i decays only into first- (second-) generation lep- 
ton. Here, we use MED (M2) model for first- (second-) generation case and 
take to 3/2 = 300 GeV, 600 GeV, and 1.2 TeV, with r 3/2 = 2.0 x 10 26 sec, 
1.1 x 10 26 sec, and 8.6 x 10 25 sec (9.3 x 10 25 sec, 5.8 x 10 25 sec, and 
5.0 x 10 25 sec), which are the best-fit lifetime with PAMELA data, respec- 
tively. Here, in order to determine the best-fit lifetime, we calculate \ 2 by 
the use of PAMELA data. (In our \ 2 analysis, since the positron fraction in 
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Fig. 2. Total flux: <E> e+ + <3? e _ with MED (M2) model for the case that the i/> M dom- 
inantly decays to the first- (second-) generation lepton in (a) ((b)). Dot-dashed line is 
the background flux. Here, we use the same parameters in (a) and (b) of Fig. 1, respec- 
tively. In (b), we also plot the flux with rri3/2 = 2 TeV and T3/2 = 4.6 X 10 25 sec and 
PPB-BETS data 24 . 



the low energy region is sensitive to the background fluxes, we only use the 
data points with E > 15 GeV.) From the figure, we see that the positron 
fraction well agrees with the PAMELA data for 7773/2 > 100 GeV irrespec- 
tive of the gravitino mass if t 3 / 2 is properly chosen. (Simultaneously, the 
energetic 7-ray flux is also enhanced, which can be an explanation of the 
7-ray excess observed by the EGRET. 11 ) 

Next, we move on to the total flux: $ e + + <£ e - . The numerical results are 
shown in Fig. 2. Here, we use the best-fit lifetime with the PAMELA data 
for each 7713/2, namely the same value in (a) and (b) of Fig. 1, respectively. 
From the figure, we see that the observed anomalous structure is well re- 
produced in the both cases. Especially, the result is a good agreement with 
the observation when m 3 / 2 = 1.2 TeV (2 TeV) for the case that the final 
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Fig. 3. Synchrotron radiation fluxes at v = 22 GHz as functions of gravitino mass for 
angle 9 = 5°, 10°, 15°, and 20°. The final-state lepton in the i/>^ decay is in the first 
generation. Here, we take r 3 / 2 = 5 X 10 26 sec, and show the cases of B =1, 3, 10 fiG 
(from the bottom to the top) for each figure. 



state lepton is the first- (second-) generation. We also note that, in the total 
flux, the numerical results does not change drastically by the choice of the 
background. This is because the signal from the gravitino is larger than (or 
at least comparable to) the background. 

Finally, let us discuss the synchrotron radiation flux. The numerical 
results are shown in Fig. 3. In this figure, we consider the case that the ^ 
mainly decays to first generation lepton and plot for v = 22 GHz as the 
function of m 3 / 2 , taking t 3 / 2 = 5 x 10 26 sec. The angle is set as <f> = f , and 
9 = 5°, 10°, 15°, and 20°, and wc take B = 1, 3, and 10 fj,G. For the W 
DM case, it can be seen that the synchrotron radiation flux is of the order 
of ~ 1 kjy/str or smaller. As we mentioned, since the the existence of the 
exotic radiation flux of this size is controversial, it is difficult to confirm or 
exclude the present scenario without better understandings of the sources 
of Galactic foreground emission. 



6. Conclusions 

In this paper, we have studied the cosmic-ray fluxes from the ^-DM de- 
cay in RPV, motivated by the recent observations by PAMELA and ATIC. 
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Assuming that the ip^ is the dominant component of DM, we calculate the 
cosmic-ray e^, and found that the both anomalies can be well explained 
when T 3 / 2 ~ O(10 26 sec). In particular, we saw that the ATIC anomaly 
indicates TO3/2 ^1 — 2 TcV in this scenario. We also calculate the syn- 
chrotron radiation induced by the cosmic-ray e from the Galactic center 
region with the lifetime to explain PAMELA and ATIC anomalies. Then, 
we obtained the result that the synchrotron radiation flux is 0(1 kJy/str) 
or smaller, which does not exclude our scenario by the observation of the 
Galactic foreground emission. 
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